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We study the neutrino mass hierarchy at the magnetized Iron CALorimeter (ICAL) detector at India-based 
Neutrino Observatory with atmospheric neutrino events generated by the Monte Carlo event generator Nuance. 
We judicially choose the observables so that the possible systematic uncertainties can be reduced. The resolu- 
tion as a function of both energy and zenith angle simultaneously is obtained for neutrinos and anti-neutrinos 
separately from thousand years un-oscillated atmospheric neutrino events at ICAL to migrate number of events 
from neutrino energy and zenith angle bins to muon energy and zenith angle bins. The resonance ranges in terms 
of directly measurable quantities like muon energy and zenith angle are found using this resolution function 
at different input values of 613. Then, the marginalized % 2 s are studied for different input values of 9b with 
its resonance ranges taking input data in muon energy and zenith angle bins. Finally, we find that the mass 
hierarchy can be explored up to a lower value of ©13 ~ 5° with confidence level > 95% in this set up. 

PACS numbers: 14.60.Pq 



I. INTRODUCTION 

The recent experiments reveal that neutrinos have mass and 
they oscillate $M from one flavor to another flavor as they 
travel. In the standard framework of oscillation scenario with 
three active neutrinos, the mass eigenstates |v,) and the flavor 
eigen states |v a ) are related by a mixing matrix U 



= Uai\Vj 



(1) 



The flavor mixing matrix in vacuum U (PDG representation 
fill) is 
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In this oscillation picture, there are six parameters: two 
mass squared differences Am?- (Amf 2 , Am 2 2 ), three mixing an- 
gles (012,823,813), and a CP-violating phase 8. At present, 
information is available on two neutrino mass-squared differ- 
ences and two mixing angles: from atmospheric neutrinos one 
gets the best-fit values with 3a error |Am 2 2 | ~ 2.50^ 52 )0 x 

1(T 3 eV 2 , 9 23 ~ 45.0°+i 4 ir; while solar neutrinos tell us 
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The present experiments are sensitive to I Am 2 2 1 only 
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and the both sign of Am^ fit the data equally well. 

The neutrino mass hierarchy, whether normal (to 2 > to 2 ), or 
inverted (to 2 > to 2 ), is of great theoretical interest. For exam- 
ple, the grand unified theory favors normal hierarchy. It can 
be qualitatively understood from the fact that it relates leptons 
to quarks and quark mass hierarchies are normal. The inverted 
mass hierarchy which is unquarklike, would probably indicate 
a global leptonic symmetry ||4j]. 
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For nonzero value of 813, the hierarchy can be determined 
from the matter effect on neutrino oscillation. The contri- 
butions of this effect to the effective Lagrangian during the 
propagation through matter are opposite for neutrino and anti- 
neutrino, which depends mainly on the value of 813 and the 
sign of Ato 2 2 . So the total number of observed events is ex- 
pected to differ for normal hierarchy (NET) and inverted hi- 
erarchy (IH). The hierarchy can be more prominently distin- 
guished if the experiment is able to count neutrinos and anti- 
neutrinos separately. This needs charge identification of the 
produced leptons. 

In case of 813 = 0, the mass hierarchy can also be observ- 
able in principle even in case of vacuum oscillation due to 
nonzero value of Ato 21 |5|]. 

Presently over the world, there are many ongo- 
ing and planned experiments: MINOS|6], T2K0], 
ICARUS [8], NO\Aji], D-CHOOZlfToll, UNOlfll, 
SKffllfTfT OPERAfl4ir Hyper-Klfll, and many others. 
It is notable that all these experiments are planned in the 
north hemisphere of the Earth. Among them MINOS is the 
only magnetized one and has a good charge identification 
capability. A large magnetized ICAL detector is under 
strong consideration for the proposed India-based Neutrino 
Observatory (INO)[15] near the equator. Since it has high 
charge identification capability (<; 95%) 11611 . it is expected 
to have reasonably better capability to determine the neutrino 
mass hierarchy. 

It is also notable here that the CERN-INO baseline happens 
to be close to the so-called 'magic' baseline fT7ilT8Tl for which 
the oscillation probabilities are relatively insensitive to the yet 
unconstrained CP phase. This permits such an experiment to 
make precise measurements of the masses and their mixing 
avoiding the degeneracy issues [19j| which are associated with 
other baselines. 

The mass hierarchy with atmos phe ric neutrinos at the mag- 
netized ICAL has been studied in H20l [2lt [22h . In I20H2II1 , the 
number of events are calculated in energy and zenith angle 
bins. Then the chi-square is calculated including the effect of 
possible uncertainties over a large range of zenith angle and 
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energy. It was found in l20ll that a precision of ~ 5% in neu- 
trino energy and 5° in neutrino direction reconstruction are re- 
quired to distinguish the hierarchy at 2a level with 200 events 
for sin 2 20i3 =0.1 and roughly an order of magnitude larger 
event numbers are required in case of resolutions of 15% for 
the neutrino energy and 15° for the neutrino direction. How- 
ever, in [21], the authors are able to relax the resolution width 
to 10 — 15% for energy and 10° for direction reconstruction in 
distinguishing the type of hierarchy at 95% CL or better for 
> 0.05 with 1 Mton-yr exposure of ICAL. 
In l22ll . an asymmetry parameter is considered as 
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A N (L/E) = -(L/E)- 
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where, u, (if), d, (d) are the number of up going and down 
going events for neutrino (anti-neutrino), respectively for a 
baseline L and an energy E. For down going events, the 'mir- 
ror' L is considered, which is exactly equal to that value when 
neutrino comes from the opposite direction. Here the asym- 
metry is calculated integrating over E > 4 GeV. Here a sta- 
tistically significant region with maximum number of events 
corresponds to the range 500 < L/E < 3000. However, the 
results depends crucially on the L/E resolutions. 

In M23I1 . the neutrino and the anti-neutrino events for direct 
as well as inverted hierarchy are considered for the E range 
5-10 GeV and the L range 6000 - 9500 km. 

In our work, we find the resonance ranges in terms of di- 
rectly measurable parameters 0™ nlth and E^ for a given set of 
oscillation parameters, where the matter effect contributes sig- 
nificantly in distinguishing the hierarchy. We have also stud- 
ied the possible observables and find the suitable ones for dis- 
crimination of the hierarchy. Finally, we make a marginal- 
ized % 2 study over a range of oscillation parameters and find 
how low 9i3 can be explored in discriminating the hierarchy 
atlNO. 



the plates. This type of detector has a very good time (~ 1 ns) 
as well as spatial resolution. 



III. ATMOSPHERIC NEUTRINO FLUX AND EVENTS 

The atmospheric neutrinos are produced from the interac- 
tions of the cosmic rays with the Earth's atmosphere. The 
knowledge of primary spectrum of the cosmic rays has been 
improved from the observations by BESSJH] and AMSll25ll. 
However large regions of parameter space have been unex- 
plored and they are interpolated or extrapolated from the mea- 
sured flux. The difficulties and uncertainties in calculation of 
the neutrino flux depends on the neutrino energy. The low en- 
ergy fluxes have been known quite well. The cosmic ray fluxes 
(< 10 GeV) are modulated by solar activity and geomagnetic 
field through a rigidity (momentum/charge) cutoff. At higher 
neutrino energy (> 100 GeV), solar activity and rigidity cutoff 
are irrelevant ll26ll . There is 10% agreement among the calcu- 
lations for neutrino energy below 10 GeV because different 
hadronic interaction models are used in the calculations and 
because the uncertainty in cosmic ray flux measurement is 5% 
for cosmic ray energy below 100 GeV j26ll . In our simulation, 
we have used a typical Honda flux calculated in 3-dimensional 
scheme ll26ll . 

The interactions of neutrinos with the detector material 
are simulated using the Monte Carlo model Nuance (version- 
3) 112711 . Here the charged and neutral current interactions are 
considered for (quasi-)elastic reactions, resonance processes, 
coherent and diffractive, and deep inelastic scattering pro- 
cesses. 

In fig. Q] we show the variation of total Charged Current 
(CC) cross section with neutrino energy for both neutrino and 
anti-neutrino. In the third plot of this figure, we also show 
how rapidly the neutrino flux changes with energy. 



II. THE INO DETECTOR 



IV. NEUTRINO OSCILLATION THROUGH THE EARTH 
MATTER 



The detector is a rectangular parallelepiped magnetized 
Iron CALorimeter detector! 1 511 ■ Here the simulation has been 
carried out for 50 kTon mass with size 48 mx 16mx 12 m. 
It consists of a stack of 140 horizontal layers of 6 cm thick 
iron slabs interleaved with 2.5 cm gap for the active detector 
elements. For the sake of illustration, we define a rectangu- 
lar coordinate frame with origin at the center of the detector, 
x(y)-axis along the longer (shorter) lateral direction, and z- 
axis along the vertical direction. A magnetic field of strength 
w 1 Tesla will be applied along +y-direction. The resistive 
plate chamber (RPC) appears to be the best option for the ac- 
tive part of the detector. It is a gaseous detector consisting 
of two parallel electrodes made of 2 mm thick 2 m x 2 m 
glass/Bakelite plates with graphite paint on the out sides and 
separated by a gap of 2 mm. When a charge particle passes 
through this active part, it gives a transient and a very local- 
ized electric discharge in the gases. The readout of the RPCs 
are the 2 cm width Cu strips placed on the external sides of 



The effective Hamiltonian which describes the time evolu- 
tion of neutrinos in matter can be expressed in flavor basis as 



H = — 
2E 



U 
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where, A = G\j2N e 2E is the effective potential of v e with 
electrons, U is the flavor mixing matrix in vacuum (eq. |2), G 
is the Fermi constant, N e is the electron density of the medium, 
and E is the neutrino energy. 

For a time t, the evolution of neutrino states is given by 
v(r) = S(t )v(0) with S(f ) = e~ iHt for constant matter density. 
The corresponding effective Hamiltonian for anti-neutrinos is 
obtained by U — ► U* and A — * —A. To understand the analyti- 
cal solution one may adopt the so called "one mass scale dom- 
inance" (OMSD) frame work: (A/n^l « \m\-m 2 X7 \ JMIHl. 
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FIG. 1: The variation of total CC iron cross section (scaling by a factor of ^) with neutrino energy for both neutrino (left) and 
anti-neutrino (right) and the variation of atmospheric flux (Honda 3D) with energy (lower) for a fixed zenith angle and azimuthal 
angle. 



With this OMSD approximation, the survival probability of 
can expressed as 



P™ = 1 - cos 2 e™ 3 sin 2 29 23 
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(5) 



The mass squared difference (Am^)™ and mixing angle 
sin 2 20|" in matter are related to their vacuum values by 



:(6) 



(Aml l ) m = A/((Am 2 1 )cos20i3-A) 2 + ((Am 2 1 )sin20i3) 
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From eq. |5]and|6]it is seen that a resonance in P™ will occur 
for neutrinos (anti-neutrinos) with NH (IH) when 

sin 2 26f 3 1 or, A = Amfj cos 20 B . (7) 



Then resonance energy can be expressed as 



1 



2 x 0.76 x 10- 4 F, 
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Fig. |2] shows how the resonance energy and the average 
density of the Earth changes with zenith angle. 



V. ROLE OF HIERARCHY IN OSCILLATION 
PROBABILITY AT DIFFERENT E v AND COS 0; ENITH 



The significant difference in survival probabilities between 

NH and IH arises due to matter resonance. For a neutrino 

energy (E v ), the corresponding baseline (L v ) or zenith angle 
(9 zenith) where the 

resonance occurs, can be understood from 
fig. |2] for both cases of v^, and v^. There are two zones of 
energy: 

1 . low energy zone (E as 4 GeV) for the events through the 
core of the Earth, and 

2. high energy zone (E w 6 — 8 GeV) for the events 
through the mantle of the Earth. 

Here we have studied in detail the difference in survival prob- 
abilities at different E v and cos 0™ nlth using full three flavor 
oscillation formula with the Earth density profile of PREM 
model[30]. The dependence on the values of the oscillation 
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FIG. 2: The variation of average density (y2-axis) and the corresponding resonance energy (y1-axis) with zenith angle (x1-axis) 
and baseline (x2-axis) at 13 resonance assuming one mass scale dominance approximation. We set 613 = 10°. 



parameters are also studied in detail. The observations are the 
following. 

The survival probabilities for both NH and IH without mat- 
ter effect produce a sinusoidal behavior in L for a fixed E. In 
case of non-zero 813 with NH, it is seen that when Vu with 
energy w 2 — 6 GeV passes through the core of the Earth 
(density w 12gm/cc), a large depletion in survival probabil- 
ity P(Vjj <-> Vf,) with respect to IH arises due to the increase 
of PiVft <-» V e ) (see first row of fig. [3}. In case of Vu, it hap- 
pens for IH (see second row of fig. [3j. The important point 
is that the atmospheric neutrino flux is sufficiently high at this 
energy range. 

Again, there is also a large depletion in survival probabil- 
ity of neutrino (anti-neutrino) for NH (IH) with respect to 

IH (NH) at the high energy (w 5 - 10 GeV) for -0.75 > 
cos0 zenith > _ 40 as seen in fig g] and [5] 

This pattern remains almost same over the present allowed 
range of oscillation parameters. 

Though all the above effects diminish rapidly as 613 goes 
to zero, there remains a finite difference in the survival prob- 
abilities for NH and IH at 813 =0 due to nonzero value of 



VI. POSSIBLE SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties may enter into the analysis of 
data in three steps: i) flux estimation, ii) neutrino interaction, 
and iii) event reconstruction. 

We can divide each of them into two categories: I) overall 
uncertainties (which are flat with respect to energy and zenith 
angle), and II) tilt uncertainties (which are function of energy 
and/or zenith angle). 



The uncertainty in flux of category II may arise due to the 
tilt in its shape with energy and zenith angle. This can be 
expressed in the following fashion for the case of energy: 



E \ 

<t> 6/: (£) =*„(£)( —J «<*>(,(£) 



l+8 £ log 10 — 
Eq 



(9) 



This arises due to the uncertainty in spectral indices. The un- 
certainty 5f =5% and Eq = 2 GeV is considered in analogy 

with SB. 

Similarly the flux uncertainty as a function of zenith angle 
can be expressed as 



4> 8 .(cos8 z ) w ^(cosSj) [1 +8 z cos6,] 



(10) 



However, this uncertainty is less (« 2%) H33L 13411 than the 
energy dependent uncertainty. 

The experimental systematic uncertainties may come 
through reconstruction of events. Till now it is not estimated 
by the simulation of ICAL at INO. However the fact from sim- 
ulation is that the wrong charge identification possibility is al- 
most zero for the considered energy and zenith angle ranges. 
Since the detector is symmetric to up and down going events, 
it is expected that the up/down ratio cancels most of the ex- 
perimental uncertainties. The quantitative details of the un- 
certainties of category I (except those arising from the event 
reconstruction) can be found in l35ll . 



VII. MIGRATION FROM NEUTRINO TO MUON 
ENERGY AND ZENITH ANGLE 

In the CC processes, a lepton of same leptonic family of 
neutrino is produced with addition of some hadrons. Among 
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FIG. 3: The variation of the (first row) and v M (second row) survival probability with cos0 zen ith with a fixed energy E v = 3.7 GeV 
for On = 4° (first column), 7° (second column), and 10° (third column). The other oscillation parameters are set at their best-fit 
values. 



the leptons, only muon produces a very clean track in ICAL 
detector and the hadrons produce a few hits around the vertex 
of the interaction. Our analysis considers only muons, and no 
hadrons for simplicity. 

For a given neutrino energy, the scattering angle and the 
energy of the muon are related to each other. Again, the neu- 
trino direction is specified by two angles: the polar angle and 
the azimuthal angle. Therefore, the actual resolution function 
should be a function of energy and these two angles. 

We are only interested in zenith angle which is a compli- 
cated function of above two angles. We reduce the above 
complexity of two angles by constructing a two dimensional 
resolution function R(xe, x Z enith) for a given neutrino energy 
with 

(Ey —E^j) _ / zenith n zenith\ /iix 

XE — ~ , ^zenith — (,t) v — H ). (11) 

Ev 

To find the resolution function R(xE,x zftn i t h), we generate one 
thousand years un-oscillated atmospheric neutrino data for 
ICAL. Here, the two angles automatically come in proper way 
in zenith angle resolution and it should work well in analysis 
of the atmospheric data. We divide the whole data set into 
17 neutrino energy bins (for E =0.8-18 GeV) in logarithmic 
scale. It is also notable that the zenith angle resolution slightly 
depends on the value of neutrino zenith angle due to the spher- 
ical shape of the Earth. To account for this dependence, we 
divide the whole zenith angle range (cos Q z = — 1 to +1) into 
10 bins for every energy bin. In fig. [6]we show typical resolu- 



tion functions at low (« 1 GeV) and high (w 15 GeV) energy 
for both neutrinos and anti-neutrinos. Here, we see, the zenith 
angle resolution improves and energy resolution worsens with 
increase of energy. The zenith angle resolution is significantly 
better for v A ,s than v^s. 

From the distribution of resolution function (fig. [6j, it is 
clear that the probability is high near zero of Jc zen ith an d it 
rapidly falls as we go far from it. Therefore we use a vary- 
ing bin size in our analysis for both variables, small near zero 
and gradually bigger at far values. 

How we obtained the number of events in muon energy and 
zenith angle bins from atmospheric neutrino flux using the 
above resolution function is described in the following steps. 

1. For a given neutrino energy and zenith angle bin, we 
obtain a distribution of events in a plane with axes 
x E = (E v -E„)/E v andxzenith = (6 v enith - 6^) froma 
large number of events generated by Nuance as shown 
in fig. 6. From this distribution we find the probabil- 
ity of getting a muon at a given energy and zenith angle 
bin dividing the number of events of this bin by the total 
number of events generated for the particular neutrino 
energy and zenith angle. 

2. The number of events for a given neutrino energy and 
zenith angle bin is calculated by multiplying the ex- 
posed flux with the cross section of this energy. (The 
cross section is obtained from Nuance as shown in Fig. 

ni) 
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FIG. 4: The same plots of fig. fj but with E v = 7.5 GeV. 



3. To obtain how many events at a particular muon en- 
ergy and zenith angle bin may come from the number 
of events of a given neutrino energy and zenith angle 
bin, we multiply this previously calculated number of 
events for the given neutrino energy with the probabil- 
ity at this muon energy and angle bin obtained in step 
1. 

4. The total contribution at a particular muon energy and 
zenith angle bin is obtained by adding the contribution 
from all neutrino energy and zenith angle bins. Finally 
we obtained the total muon events at different energy 
and zenith angle bins in the considered ranges. 



VIII. CHOICES OF OBSERVABLES AND ITS VARIABLES 

The effect of matter on the oscillation probability is a com- 
plicated function of E, L and the density of medium (p). From 
the study of survival probability P(v^ «-> V^) and /"(v^ <-> V^) 
in section[V] it is clear that the difference between NH and IH 
is more significant for particular ranges of E and L rather than 
the range of the combination L/E. In case of the combination 
L/E, there remains some L(E) values for some values of E(L) 
where the difference between NH and IH is insignificant. For 
these values of L and E, it increases the magnitude of statis- 
tical errors and eventually kills the significance which comes 
from the interesting regions of L and E. For this reason, we 
need to consider only important ranges of the baseline L, for 
each range of energy Ej. 



We may consider the following observables to determine 
the hierarchy: i/dt, Ui/di, or any combina- 

tions of them, where, m, , dt («,•, di) are the up going and down 
going neutrino (anti-neutrino) events in the interesting L\ — E\ 
regions. For down going events we consider the 'mirror' L 
which is exactly equal to that value when the neutrino comes 
from the opposite direction. 

In case of observables with the combination of neutrinos 
and anti-neutrinos, we will encounter the uncertainty of the 
ratio V/j/Vu which is ~ 1 at low energy and it increases with 
energy 12611311, 13211 . Since the matter effect and the cross sec- 
tions for neutrinos and anti-neutrinos are dependent on E, the 
uncertainty in the V^/v^ ratio may affect significantly in de- 
termination of hierarchy for the combined observables. 

Firstly, we choose the observables as a ratio of u and d to 
cancel all the overall uncertainties (discussed in section [VTb . 
Secondly, since the magnetized ICAL detector has the capa- 
bility to distinguish /j~ and we consider the ratios sepa- 
rately for v^s and v^s: 

A = u/d, and A = u/d. (12) 

In fig. |2j it clearly shows that there is a sudden fall in res- 
onance energy due to a sudden rise in average density of the 
Earth. We choose two ranges of energy: one for core (vertical 
events) and another for mantle (slant events). The difference 
in survival probability between NH and IH can also be seen in 
figs. [3] and|4]at these two ranges of zenith angle for two typ- 
ical resonance energy of 3.7 GeV and 7.5 GeV, respectively. 
The resonance range of energy for a typical value of the zenith 
angle is also shown in fig. [5] 
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FIG. 5: The variation of the v p (first row) and (second row) survival probability with energy for 9 13 
column), and 10° (third column) with a typical zenith angle ze nith = 130°. 



= 4° (first column), 7° (second 



From the discussion of section Hill M and TVIII we observed 
the following : 

1 . The zenith angle resolution width is large with a peak 
around zero (fig. |6). 

2. From the distribution of the energy resolution, it is clear 
that a muon of energy Ef, can come from any neutrino 
having energy E v > Ef, with almost equal probability for 
E v ^ 3 GeV where the deep inelastic events dominate. 
Again, as the flux falls very rapidly with energy (see 
fig. [TJ and E v > £„ (this is obvious for E v > 1 GeV for 
targets at rest), the lower value of the energy range will 
play a crucial role in the results. 

The bonus point of these two ranges of zenith angle (verti- 
cal and slant) is that the vertical/horizontal uncertainty in flux 
is minimized and becomes negligible. Among the systematic 
uncertainties discussed in section [VT1 we remain with the en- 
ergy dependent flux uncertainty. Here we will not consider 
any uncertainty coming from reconstruction of the events. 
This will be taken care of in near future in GEANT-based 
studies. 



IX. x 2 ANALYSIS 



It should be noted here that the resonance width falls very 
rapidly with 613. 

From the earlier works 1I20I l2ll I22I1 it is found that as the 
energy and the angular resolution width increase the discrim- 
ination ability of hierarchy becomes worsened (see fig. 10 of 
lHj]]). In our actual simulation it is found that the width of 
these resolutions are large (see fig. 6). 

We have tried to tackle the problem in the following way. 
We do not divide the selected ranges. We calculate the total 
up-going (uj) and the total down-going (di) events and then 
their ratio (ui/df) for each selected range. Finally we perform 
the chi-square study with these ui/di ratios of all ranges (see 
table|B- 

For a set of input values of the oscillation parameters, we 
generate the neutrino events for 1 MTon.year exposure of 
ICAL detector using the Nuance (see fig. |7). Next we cal- 
culate the total iii and the total c/, for each considered range of 



E u: — cos 9 



zenith 



(see tableU) and then their ratio iij/di. We call 
it the "experimental data". The statistical error is estimated by 
the following formula 



di \ di 



(13) 



We adopt the pull method for chi-square analysis. The the- 
oretical data in chi-square is calculated with the following nu- 
merical method. 



Here we are interested to find the hierarchy discrimination 
ability of ICAL@INO at low 813 using the resonance ranges. 



1. Considering a set of the oscillation parameters and a 
value of the pull variable for a considered systematic 
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FIG. 6: The distribution of probability of the resolution function for neutrino (left) and anti-neutrino (right) at low energy 1 GeV) 
(top) and at high energy (fs 15 GeV) (bottom). 
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FIG. 7: The typical distribution of the number of events with cose" ^ 
set Amj 2 = -2.5 x l(T 3 eV 2 , 6 23 = 45° and 6i 3 = 7.5°. 



for fixed E p (left) and with E^ for fixed cose? ., (right). We 



uncertainty (using eq. [9j, we obtain the neutrino flux 
in 200 bins of E v and 300 bins of cos9f nith . At low 
energy region the number of oscillation swings with E 
(L) for a fixed L (E) are very large (see fig. [3]and|5]l. 
Again the difference of the oscillation probabilities be- 
tween NH and IH changes very rapidly for some ranges 
of L and E. To obtain the accurate oscillation pattern 



one needs finer binning of L and E. We have checked 
the oscillation plots varying the number of bins for E v 
and cos9f nith and finally fix 200 bins for E v and 300 
bins for cos 

Qzenith We have use( j 

same number of bins 
to obtain both experimental and theoretical data. The 
number of events from neutrino to muon energy and 
zenith angle bins are obtained using the resolution func- 
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tion and the cross section as discussed earlier. 

2. Next we calculate the total m,- and the total of, separately 
for each selected range of E u ,-cos 

gzenith (see table 

Then the ratio ui/di is determined for each range. We 
call it the "theoretical data". It should be noted that the 
number of bins has been changed in chi-square from 
that of the flux. (The reason is described above.) 

3. We calculate the standard chi-square from these "the- 
oretical data", the above "experimental data" and the 
"statistical error" (eq. fl~3T > for all ranges. Next the min- 
imization of this chi-sqaure is done with respect to the 
pull variable to find its best-fit value keeping all other 
oscillation parameters fixed. Finally the "theoretical 
data" is again calculated using the best-fit value of the 
pull variable. 

In the last the % 2 is calculated from the final "theoretical 
data" and the above "experimental data". The whole process 
is done to find the % 2 and %j, for v^s and v^s, respectively. At 
the end, we find the total % 2 — % 2 + for the considered set 
of oscillation parameters. 

As discussed in section IVTl we consider only neutrino flux 
uncertainty due to its tilt with energy ||37ll . 

X. RESONANCE RANGES IN TERMS OF £ p AND 6J enith 

The optimization of the ranges by analytical calculation is 
almost impossible for the complexity in the behavior of sur- 
vival probability, the smearing of resolution, and the rapid fall 
of flux with increase of energy. 

We calculate the % 2 varying the energy and zenith angle 
ranges for both NH and IH for a given input of 813 and type 
of hierarchy keeping all other oscillation parameters at their 
best-fit values. Then, we select the ranges where the % 2 differ- 
ence for NH and IH is the largest. This is done for the input 
values of 813 = 4°, 5°, 7°, 10° and tabulated in tableU Due to 
the effect of energy resolution function and the nature of flux, 
the higher end of the energy range can not be precisely deter- 
mined. We see, the ranges squeeze rapidly with the decrease 
of 813. It should be noted that the ranges determined here has 
a significant weight of the atmospheric neutrino flux and it 
is not independent of flux since it is not from the migration 
of neutrinos to muons only. However, there is no significant 
change in chi-square if the ranges are changed within few per- 
cent of their boundary values. 

XI. MARGINALIZATION AND RESULTS 

The neutrino oscillation parameters can not be measured at 
infinite accuracies. To find at what statistical significance the 
wrong hierarchy can be disfavored, we calculate the % 2 for 
both true and false hierarchy varying the oscillation parame- 
ters in the following ranges (unless we specify) : 

1. \A m y : 2.3 x l(T 3 -2.8 x l(T 3 eV 2 , 



2. 893 : 39° -51°, and 

3. 8b : 0°- 12°. 

We set other oscillation parameters at their best-fit values and 
choose CP-violating phase 8 = 0. For the data set of each in- 
put values of 813, we used the corresponding resonance range 
ofE /i -cos9£ enith . 

The marginalized % 2 for both true and false hierarchy as a 
function of 813 is shown in fig. [8] for the input of IH and dif- 
ferent values of 613 (5°, 7° and 10°). Here the marginalization 
is done over other two oscillation parameters (Am^l and 823. 
The similar plot is also shown for the input of NH in fig. [9] but 
with Ep — cos 8 zemth ranges as optimized for IH. Finally we 
show the discrimination sensitivity A% 2 = % 2 (false) — % 2 (true) 
marginalized over all three oscillation parameters |Am| 2 |, 823 
and 813 as function of input values of 813 in fig. |T0j 

The changes reflected in % 2 with the input of 813 and the 
type of hierarchy can be understood from the following facts: 

1 . Normally the number of events is larger for v^, than 
due to larger v^, cross section. So the Vu contribution 
to the difference in chi-square due to hierarchy is larger 
than Vfj. Since v^s are suppressed in NH, the input with 
IH (larger statistics) will give a larger difference in chi- 
square and NH can be more easily disfavored than IH. 
(See fig. [8]and|9]for the input value of 813 = 10°.) How- 
ever, once any one type of the hierarchies is disfavored, 
there remains the other type only. So it is essential to 
disfavor any one of them. 

2. The resolution is better for anti-neutrinos than neutri- 
nos. Again as 813 decreases, the resonance ranges 
squeeze very rapidly. This leads to a competition be- 
tween resolution and statistics resulting a comparable 
contribution of and v A , to chi-square difference for 
013 ^ 7° (compare fig. [8] and [9] for the input value of 
013 = 7°). 

3. Though there is an appreciable difference in survival 
probabilities between NH and IH below 813 value of 
5°, the tilt uncertainty with energy kills the difference 
in chi-squares. 



XII. DISCUSSION 

Here we will discuss some points for the estimation of de- 
tector efficiency and the errors that may occur in event recon- 
struction of the experimental data. 

1. One can estimate the efficiency of the ICAL detector 
in the following procedure. The fiducial volume of the 
detector can be estimated so that an event with vertex 
inside this volume can be reconstructed. To reconstruct 
an event one may need the minimum number of hits 
<; 9. Therefore, one can estimate the efficiency just ne- 
glecting 10 layers from top and bottom and 50 cms from 
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Resonance ranges in terms of muon energy and zenith angle 


G13 


Core 


Mantle 


cos6 zenith 

/J- [1+ 


(GeV) 


cos6 zenith 


E^ (GeV) 
[T n+ 


10° 

7° 
4° 


-1.0 to -0.57 -1.0 to -0.72 
-1.0 to -0.63 -1.0 to -0.79 
-1.0 to -0.50 -1.0 to -0.83 


0.79 to 9.40 0.93 to 7.40 
1.39 to 10.18 2.07 to 4.59 
1.63 to 8.68 2.07 to 6.83 


-0.57 to -0.41 -0.68 to -0.46 
-0.61 to -0.43 -0.70 to -0.46 
-0.48 to -0.30 -0.72 to -0.46 


3.33 to 8.01 5.38 to 8.68 
4.24 to 6.83 5.38 to 11.94 
4.59 to 6.31 7.40 to 10.18 



TABLE I: The resonance ranges of cos 9£ enith - for different input values of 613. The type of input hierarchy is inverted. 




2 4 6 8 10 12 2 4 6 8 10 12 14 2 4 6 8 10 12 14 



FIG. 8: The variation of marginalized x 2 with e B for NH and IH for the input of IH and the values of e 13 =5° (left), 7° (middle) 10° 
(right). The marginalization is done over [Am^ 2 | = (2.3-2.7) x 10~ 3 eV 2 and 9 2 3 = 39° -51°. We consider a uncertainty in tilt factor 
8 £ = 5% at E v =2 GeV. The other oscillation parameters are set at their best-fit values (discussed in the text). 



lateral sides. For this into consideration, the efficiency 
may be roughly 

^ 120 (layers) x 4700 cm x 1500 cm _ ^ 
~ 140 (layers) x 4800 cm x 1600 cm ~ 

2. For the maximum difference between NH and IH, we 
find the lowest value of » 1 GeV at cos0£ enith = -1 
and the highest value of £ A , w 15 GeV at cos0^ emth « 
—0.3. For these cases, we will get number of hits ^ 12, 
and there will be a very high efficiency of filtration and 
reconstruction of these events in actual experiment. 

3. For the muons below the energy of 15 GeV, the charge 
identification capability of ICAL@INO is > 95% with 
a magnetic field of 1 Tesla ll 1511 . Here, we consider 
only the muons whose reconstruction efficiency is very 
high and the energy and angle resolutions are within 
3-6% for the considered ranges. Moreover, the wrong 
charge identification possibility is very negligible [15]. 
So one may easily expect that the result esti- 
mated in this paper will not change appreciably af- 
ter GEANT[36] simulation. 



trino events generated by the Monte Carlo event generator 
Nuance(version-3). We have done the analysis in muon en- 
ergy and zenith angle bins. We have adopted a numerical 
technique for the migration of number of neutrino events from 
neutrino energy and zenith angle bins to muon energy and 
zenith angle bins to find the "theoretical data" in chi-square 
analysis. The so called "experimental data" are obtained in 
muon energy and zenith angle bins from the Nuance simu- 
lation. The resonance ranges in terms of muon energy and 
zenith angle are found for different values of 613. Then we 
find suitable variables and make a marginalized chi-square 
analysis using the pull method for both true and false hier- 
archy. The choice of variables and the pull method minimize 
the possible systematic uncertainties in the analysis. Finally 
we find the sensitivity as the difference A% 2 = % 2 (false) — 
% 2 (true) as a function of On as shown in fig. |T0j It is found 
that the neutrino mass hierarchy with atmospheric neutrinos 
can be probed up to a low value of On « 5° at ICAL@INO by 
judicious selection of events and observables. Here we have 
also shown that > 95% CL can be achieved for discrimination 
of hierarchy at 9n 5° in a 1 Mton-year exposure of ICAL 
detector. 



XIII. CONCLUSION 



We have studied the neutrino mass hierarchy at the 
magnetized ICAL detector at INO with atmospheric neu- 
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FIG. 9: The same plots of fig. U but with input of NH and e 13 =7° (left) and 10° (right). 




FIG. 10: The variation of the discrimination sensitivity A% 2 with the input value of 6b for two marginalization ranges of lAm^l- The 
range of 623 is 39° - 51°. The type of input hierarchy is inverted. 
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